A bi-criterion sequence-dependent scheduling problem with order deliveries

The manufacturing sector faces unprecedented challenges, including intense competition, a surge in product varieties, heightened customization demands, and shorter product life cycles. These challenges underscore the critical need to optimize manufacturing systems. Among the most enduring and complex challenges within this domain is production scheduling. In practical scenarios, setup time is whenever a machine transitions from processing one product to another. Job scheduling with setup times or associated costs has garnered significant attention in both manufacturing and service environments, prompting extensive research efforts. While previous studies on customer order scheduling primarily focused on orders or jobs to be processed across multiple machines, they often overlooked the crucial factor of setup time. This study addresses a sequence-dependent bi-criterion scheduling problem, incorporating order delivery considerations. The primary objective is to minimize the linear combination of the makespan and the sum of weighted completion times of each order. To tackle this intricate challenge, we propose pertinent dominance rules and a lower bound, which are integral components of a branch-and-bound methodology employed to obtain an exact solution. Additionally, we introduce a heuristic approach tailored to the problem’s unique characteristics, along with three refined variants designed to yield high-quality approximate solutions. Subsequently, these three refined approaches serve as seeds to generate three distinct populations or chromosomes, each independently employed in a genetic algorithm to yield a robust approximate solution. Ultimately, we meticulously assess the efficacy of each proposed algorithm through comprehensive simulation trials.


INTRODUCTION
In today's manufacturing landscape, challenges such as intense competition, a proliferation of product varieties, heightened customization demands, and shortened product life cycles have become more pressing than ever.Consequently, there is greater focus on the need to optimize manufacturing systems.Among the well-established scheduling problems, the customer order scheduling problem (COSP) has always held significant importance in differential evolution algorithm for the same problem in de Athayde Prata, Rodrigues & Framinan (2021b).
To align with the reality of customer order scheduling, this study considers the classification of jobs from different orders into distinct classes.Every order consists of at least one job from each job class, necessitating a setup time whenever the machine transitions between classes.In the domain of single-machine scenarios with multiple objectives and various job classes, Liao (1993) introduced a branch and bound (B&B) algorithm to systematically investigate all viable solutions.Gupta, Ho & van der Veen (1997) were pioneers in investigating COSP with multiple classes and setup times, introducing a bi-criteria problem.They aimed to minimize both makespan and the total carrying costs of customer orders, with carrying costs computed based on the time difference between the completion times of the first and last job in each customer's order.Lin, Yin & Liu (2013) addressed a sequence-dependent scheduling problem that incorporates order delivery, presenting a binary integer program and a dynamic programming algorithm.They also introduced tabu search, iterated local search, and genetic algorithm approaches for obtaining approximate solutions.In related studies on customer order scheduling with setup times, Erel & Ghosh (2007) examined customer orders with varying quantities of products from different product families processed on a continuously available machine in any sequence.Liu (2009), Liu (2010) and Hsu & Liu (2009) explored various scenarios involving multiple jobs within a job shop setting.Recently, Li et al. (2023) expanded upon the work laid down by Gupta, Ho & van der Veen (1997) by introducing an innovative approach that incorporates four heuristics, three locally enhanced search methods, and a branch and bound algorithm.Additionally, they applied a water wave optimality algorithm, offering four distinct wavelength variants.This collaborative effort is aimed at minimizing a mixed criterion including the order ranges, order tardiness, and total job completion times.In a parallel effort, Lin et al. (2023) presented four heuristics complemented by a local search method, four theoretical simulated annealing techniques, a cloudy theoretical simulated annealing hyperheuristic, and a B&B method to solve the problem.Their collective endeavors focused on minimizing a mixed criterion including the total order ranges and total order tardiness.Meanwhile, Gupta et al. (2023) introduced six two-phase heuristics, six variants of water-wave optimization algorithms, a mixed-integer linear programming formulation, and a B&B algorithm to solve their problem.This concerted effort is directed at producing approximate solutions that effectively minimize a mixed criterion of the cumulative holding costs for all orders and makespan for all jobs.
Setup times present a significant concern for manufacturing companies.Lengthy setup times can jeopardize the on-time delivery of customer orders, leading to potential losses in customer satisfaction and incurring direct or indirect costs.Effectively incorporating and managing the setup time factor in scheduling decisions can yield substantial benefits (Kim & Bobrowski, 1997;Allahverdi, Gupta & Aldowaisan, 1999;Allahverdi et al., 2008;van Donk & van Doorne, 2016;Zhao et al., 2018;Ying et al., 2023).
Motivated by the scarcity of research that address sequence-dependent bi-criteria scheduling problems that incorporate order deliveries, we introduce a novel problem in this domain.Our aim is to pinpoint a schedule that optimally balances the makespan and the weighted completion time of designated orders through a linear combination.Expanding on these insights, this study addresses a sequence-dependent bi-criteria scheduling challenge involving order deliveries.The primary aim is to devise an optimized schedule that minimizes a mixed criterion of the overall job makespan and the aggregated weighted completion time of all orders.
The main contributions of this study are summarized as follows: (a) we address a sequence-dependent bi-criteria scheduling challenge involving order deliveries; (b) we propose three lemmas and a lower bound in a branch-and-bound algorithm for finding an optimal solution; (c) we propose a simple heuristic according to the smallest first value of d uv and their three improvements by three local searching methods and (d) we build three population-based genetic algorithms.The subsequent sections of this study are structured as follows.'Problem Statement' outlines the proposed problem formulation.In 'Heuristic methods and heuristic-based genetic algorithm' , we present three dominant properties, a proposed lower bound, a heuristic, and its three improved schemes (i.e., pairwise interchange, extraction and forward-shifted reinsertion, and extraction and backwardshifted reinsertion).Additionally, we propose three variants of the genetic algorithm (GA) based on these three improved heuristics.'Tuning the Related Parameter of the Three Heuristic-Based Genetic Algorithms' focuses on testing the appropriate parameter values in the GA.'Computational Simulations and Discussions' provides a comprehensive report of all test results.The final section offers conclusions and outlines future research directions.

PROBLEM STATEMENT
The formal framework of our proposed study is outlined as follows.Let us consider a set = O 1 ,...,O m comprising m orders designated for processing on a single machine.The machine cannot break down while performing all jobs.Each order O i consists of n i jobs, each assigned a weight w i (w i ∈ (0,1)).The total number of jobs across all orders amounts to n, which defines the set N = {J 1 ,...,J n }.
A critical assumption In this model (see examples in Lin, Yin & Liu (2013)) is that, irrespective of their respective orders, if job J v immediately follows job J u in the schedule, then there is a sequence-dependent setup time d uv > 0 required.Additionally, we define d 0J u to signify the machine status before the initial scheduling of job J u , elucidating the setup required.
Notably, in this model, we assume that the processing times for all jobs in set N are negligible.Given a complete schedule σ encompassing all jobs, we introduce the notation C Op (σ ) to represent the completion time for order O p , which is the time at which the last job of O p is finished.Simultaneously, the makespan, denoted as C max (σ ), signifies the completion time of the last job in set N.
The primary objective of this study is to determine an optimal schedule that minimizes a linear combination of the sum of makespan and the weighted completion times of all m orders.In mathematical terms, we aim to find a solution to the following optimization problem: Minimize g(σ Here, θ represents a weighting parameter that allows for the adjustment of the trade-off between makespan and weighted completion times. To simplify matters, we use the terms ''node,'' ''job,'' and ''city'' interchangeably to describe the fundamental elements of our scheduling problem.Additionally, we introduce two notations, C u (σ ) and C v (σ ) to represent the completion times of jobs u and v in schedule σ , respectively.C v σ and C u σ represent the completion times of jobs v and u in schedule σ , respectively, where σ = (π ,J u ,J v ,π c ) and σ = (π ,J v ,J u ,π c ).Here, π and π c are two distinct sequences within set N.
To enhance the efficiency of our B&B method, we establish two properties that assist in narrowing down the nodes that satisfy the conditions outlined in these properties.Before presenting the properties, we designate job J L as the final job in sequence π.
Then, g σ − g(σ ) > 0. Let t π be the completion time of the last job in π, and t p1 (t q1 ) be the completion time of the first job of O p (O q ).Note that t p1 < t π (t q1 < t π ).Moreover, for any possibly arranged job sequence of all jobs in π c , let J s be the first assigned job in π c .Note that J s = J u , J s = J v .Then, one has the following results, (i) and (ii).
w q +w p d vu − w q w q +w p d uv > 0. For these orders (not including O p and O q ), which are the last jobs of the order assigned in π c , the sum of all w Combining the results of (i) and (ii), g σ > g(σ ) follows.Property 2: As jobs J u ∈ O p , J v ∈ O q , p = q, if J u is not the last assigned job of O p , and J v is also not the last assigned job of O q , d Lu + d uv < d Lv + d vu and max x∈π c d vx < min y∈π c d uy , then σ dominates σ .Proof: It needs to be shown that g(σ ) < g σ .
For any possibly arranged job sequence of all jobs in π c , one can be shown that The last sum is over all orders on which the last job of the order is in π c .Thus, g(σ ) < g σ .Property 3: Consider a scheduled σ = (PS,US), where PS and US denote the scheduled part with k jobs and the unscheduled part with (n-k) jobs, and and J v and J [k] belong to the same order, then σ = (PS,US) can be eliminated.
In the ensuing discussion, our objective is to establish a lower bound for the schedule denoted as σ = (π ,π c ), where π represents the scheduled portion comprising k jobs, and π c signifies the set of n-k unscheduled jobs.Let t π denote the completion time of the final job, J L , in π, and i∈π w i C O i stand for the aggregate weighted completion time of the scheduled orders in π .
In light of the sequence-dependent aspect, we incorporate J L L into the set π c , resulting in the creation of a new set denoted as π c * = π c ∪{J L }, now comprising (n − k+1) jobs.To approximate the completion times in π c , we embark on the following procedure to extract the (n − k) smallest sequence-dependent setup times within π c * .These recorded values are denoted as D(1), D(2), . . ., D(n − k).
Notably, the values of D(1), . . ., D(n − k) exhibit a non-decreasing trend.Consequently, we can estimate the lower bound for each C [k+i] (σ ), with i ranging from 1 to n − k, by employing the following formula: Additionally, we maintain a record of the frequency number f i corresponding to each order in π c , where i spans from 1 to q, with 1 ≤ f 1 ≤ f 2 ≤ ... ≤ f q and q<m.Utilizing Eq. ( 1), we derive completion times for the outstanding orders in π c as C Furthermore, we arrange the weights w (1) ≤ w (2) ≤ ... ≤ w (q) in non-decreasing order within the set {w v ∈ π c ,v =1 ,...,q}.Applying the lemma established by Hardy, Littlewood & Polya (1967), we arrive at a lower bound for the estimated weighted completion times of the unscheduled q orders, represented as q i=1 w (q−i+1) C k+ i v=1 f v (σ ).Hence, the lower bound is derived as follows: , where b +q =m.

Heuristic methods and heuristic-based genetic algorithm
Given that the proposed problem with only a single criterion has been shown to be NP-hard, the problem require significant CPU time to find exact solutions as the number of jobs increases.Thus, three heuristics based on problem characteristics will be provided.In the beginning, we choose a job J u with the smallest value among job set N, that is, d 0J u = min{d 0J i , J i ∈ N } to be scheduled on the first position.Subsequently, we choose the second job J v with the smallest value among job set N \{J u }; i.e., d J u J v = min{d J u J i ,J i ∈ N \{J u }}, and we perform the same process until all jobs are determined to yield an initial schedule σ .To obtain a finer quality of solutions, following the idea of Della Croce, Narayan & Tadei (1996), we perform three improved schemes (pairwise interchange (PI), extraction and forward-shifted reinsertion (FOR) and extraction and backward-shifted reinsertion (BK)) on the initial schedule σ .They are denoted as DPI, DFOR, and DBK.The details of the proposed heuristics are given as follows: The procedures of three improved algorithms: However, the existing level of complexity is insufficient to address the intricacies of the proposed model.Conversely, metaheuristics present a notably higher level of complexity and present a significant challenge in terms of effective design and implementation, especially in the context of intelligent random search strategies (Holland, 1975;Bean, 1994;Larranaga et al., 1999;Nguyen, Mei & Zhang, 2017).
Within the realm of metaheuristics, the GA distinguishes itself as a well-established and effective method for generating high-quality approximate solutions to a diverse set of combinatorial problems.Extensive research has demonstrated the GA's ability to discover near-optimal solutions for a wide array of complex challenges (Iyer & Saxena, 2004;Essafi, Matib & Dauzere-Peres, 2008;Beasley, Bull & Martin, 1993;Fan et al., 2022).
The foundational principles of GA were rigorously established by Holland (1975).A genetic algorithm initiates with a set of feasible solutions, referred to as the population, and subsequently replaces the current population through an iterative process.This entails the necessity of a suitable encoding for the specific problem and the formulation of a fitness function that serves as a metric for assessing the quality of each encoded solution (chromosome or individual).The reproductive mechanism selects the parents and employs a crossover operator to recombine them, resulting in offspring that undergo a mutation 00: Input {d 0J i , J i ∈ N } and {d J i J j ,J i ,J j ∈ N } 01: Set σ as an empty set and l =1 02: Choose J u with the smallest value of d 0J u among {d 0J i , J i ∈ N } to be the l th position 03: Set σ = (J u ,...), delete J u from N, record as N \{ J u }, and l =l +1 04: Do while {l <= n} 05: Choose J v with the smallest value of d J u J v among {d J u J i ,J i ∈ N \ {J u }} to be the l th position 06: Set σ = (J u ,J v ,...), delete J v from N {J u } 07: Record as N \{J u ,J v }, and l = l + 1 08: End while.09: Output the final solution σ .10: Apply the methods (PI, FOR, and BK) separately to improve σ 11: Do k 1 =1 to n − 1 12: Do k 2 =k 1 +1 to n 13: Choose the k 1 position job and the k 2 position job in σ to perform the PI, FOR, and BK separately to improve σ .14: Keep the best one.15: End do 16: End do 17: Output these three refined solutions (DPI, DFOR, and DBK).
operator to introduce local alterations (Essafi, Matib & Dauzere-Peres, 2008).The genetic algorithm proceeds through the following key stages: Step 1. Representation of structure: A structure is defined as a sequence of jobs (Etiler et al., 2004) Step 2. Initial population: Drawing on established heuristics in the literature is advised for constructing the initial population (Etiler et al., 2004).This approach expedites the convergence toward the final solution.This study adopts three initial sequences: DPIGA, DFORGA, and DBKGA create their populations based on three initial sequences DPI, DFOR, and DBK by using the pairwise interchange method.
Step 3. Population size: The population is set at nsize.

Step 4. Fitness function:
The objective is to minimize the objective function g (σ ) (i.e., a linear combination of the sum of the makespan and the weighted completion times of all m orders, i.e., g The fitness function of each string (i.e., ith string in the t th generation) is calculated as and the probability of selection of the ith string is defined as This function is pivotal for parent selection and reproduction.
Step 5. Crossover: This study employs the partially matched crossover technique, with a crossover rate of 1 Step 6. Mutation: The mutation rate is set at p.
Step 7. Selection: Population sizes remain constant at nsize from one generation to the next.Offspring, excluding the best schedule with the minimum g (σ ), are generated from parent chromosomes using the roulette wheel method.

TUNING THE RELATED PARAMETER OF THREE HEURISTIC-BASED GENETIC ALGORITHMS
In this section, we conduct several experimental tests to explore the values of related parameters of three heuristic-based genetic algorithms (GAs).To determine the parameters used in three GAs as chromosome count nsize, evolution generations gsize and mutation rate p, we consider an experimental design as the number of jobs at n =12, the number of orders at m =2, and the weight of objective function at θ = 0.25.One hundred instances are tested for this case.The initial setup times d 0J u and sequence-dependent setup time d uv are generated from uniformly distributed U[1, 20], while the weights of orders w i are generated from uniformly distributed U[1, 10].The optimal solutions (Opt i ) and the near-optimal solutions (GAH i ) are produced by performing the B&B method and three heuristic-based GAs (DPIGA, DFORGA, or DBKGA).The measurement criterion is based on the average error percentage (AEP) or AEP = 100 i=1 GAH i − Opt i /Opt i /100.The parameter tuning process is based on a single-factor experimental method applied to a set of generated problem instances.
To evaluate the value of nsize, gsize was fixed initially at 50 and p was fixed at 0.02.We executed DPIGA, DFORGA, and DBKGA several times and changed the values of nsize over the interval [10,50] with an increment of four.As shown in Fig. 1 (top), as the nsize value increases, the trends of the AEPs rapidly decrease and become flat at point 50 obtained from DBKGA.Therefore, in subsequent simulation experiments, the nsize is adopted at 50.
With the previously determined nsize at 50 and fixed mutation rate at 0.02, we performed three GAs to determine the values of gsize over [10,200] with an increment of 10. Figure 1 (middle) shows that the AEP of DBKGA has a minimum point at 180.Thus, gsize is set at 180 for all three proposed GAs in later experiments.
After investigating the parameters nsize and gsize, we established nsize =50 and gsize =180.We conducted three GAs to assess the mutation rate p, ranging from 0.01 to 0.2 with a consistent increment of 0.01.As depicted in Fig. 1 (bottom), AEP demonstrates a gradual convergence, while DBKGA attains the minimal error at a mutation rate of 0.2.Accordingly, we adopt a mutation rate of 0.2 as the benchmark for all three GAs in the subsequent stages.
In a similar way, we conducted another one hundred random instances at n =240 to explore the values of three parameters of nsize, gsize, and p for three GAs (DPIGA, DFORGA, and DBKGA).Based on previous studies, to capture parameters that outperform three  1.It can be observed in Table 1 that nsize =100, gsize =1500, and p =0.2 are appropriate values for using DPIGA, DFORGA, and DBKGA.
Based on the above tests, the population size (nsize), generation size (gsize), and mutation rate (p) are set at (50, 180, 0.2) for the small n case and at (100, 1500, 0.2) for big n case in all three GAs (DPIGA, DFORGA, and DBKGA).

COMPUTATIONAL SIMULATIONS AND DISCUSSIONS
This section conducts some tested designs to assess the performances of the proposed B&B, three heuristic methods (DPI, DFOR, and DBK), and three heuristic-based GAs (DPIGA, DFORGA, and DBKGA).The solution methodologies for the proposed problem are summarized as follows.All the related parameters are designed as follows.The initial setup times d 0J u and sequence-dependent setup time d uv are generated from uniformly distributed U[1, 20], and the weights of orders w i are generated from uniformly distributed U[1, 10].The experimental designs include the weight of the objective function at θ = 0.25, 0.5, and 0.75, the number of jobs at n =12, 16, 20, and 24 for the small-job case, and n =60, 120, 180 and 240 for the large-number-of-job case.It should be emphasized that n =m*nc signifies that the total of n jobs can be apportioned among m orders, with each order containing nc jobs.The details are summarized in Table 2.The experiment investigates 24 cases for both small and large job numbers.In each case, 100 instances are generated to assess the performance of all proposed methods.In total, 2,400 instance problems are examined for both small and large job sets.The findings are summarized in Tables 3, 4, and 5.
When dealing with a case with a small number of jobs, we assess the effectiveness of the three heuristic algorithms (DPI, DFOR, and DBK) and the DPIGA, DFORGA, and DBKGA algorithms by examining their average error percentage and maximum error across various values of n, m, and θ .The AEP is computed as follows: AEP = 100 1 Hsol i − Opt i /Opt i /100, where Opt i is derived from the B&B method, while Hsol i isproduced using the six proposed approximate techniques.Additionally, we furnish the mean and maximum values of the explored nodes, along with their respective CPU times (measured in seconds), to evaluate the performance of B&B for each scenario.It is worth noting that an instance problem is considered feasible if the B&B process requires fewer than 10 8 nodes.
With regard to the performance of B&B, columns 3 and 5 of Table 3 show that B&B consumes more nodes or CPU time at smaller order numbers (m =2) compared to larger order numbers (m =3 or 4).This implies that the properties and lower bounds become more useful as the number of orders increases.The effect of θ in columns 3 and 5 of Table 3 does not show a clear trend.Additionally, the number of nodes increases significantly as the number of jobs (n) increases.This phenomenon aligns with the characteristics of the NP-hard problem.
In regard to the obtained results for the three heuristics (DPI, DFOR, and DBK) and three genetic algorithms (DPIGA, DFORGA, and DBKGA) applied to cases with small n, Table 4 and Fig. 2 present their respective average error percentages.To summarize, the mean AEP values are 0.4733, 0.5246, and 0.4551 for DPI, DFOR, and DBK and 0.3412, 0.3504, and 0.3358 for DPIGA, DFORGA, and DBKGA, respectively.The DBKGA algorithm performed better than all other approaches.Table 4 provides the mean AEP values corresponding to different values of parameters θ , m, and n.Upon inspection of Table 4, it is evident that the AEPs for the three heuristics and three GAs experience a marginal rise as n increases  from 12 to 24.In contrast, these AEPs exhibit a comparatively stable trend when subjected to changes in the values of m or θ .According to the criterion of average relative percentage deviation (RPD), RPD is calculated as follows: RPD = 100 1 (Hsol i − H min )/H min /100, where Hsol i represents the output obtained using the six proposed approximate methods, and H min is the minimum value among these six solutions.Table 5 and Fig   Regarding computational times, Fig. 4 displays violin plots that visualize the distribution of average CPU time for both the three heuristic algorithms and the three GAs under the conditions of large n.These plots clearly indicate that the overall range and interquartile ranges of average CPU times for DPIGA, DFORGA, and DBKGA are broader compared to DPI, DFOR, and DBK.Specifically, the average CPU time values are 0.0196, 0.0195, and 0.0197 for DPI, DFOR, and DBK, respectively, while they are 2.5596, 2.6028, and 2.5983 for DPIGA, DFORGA, and DBKGA, respectively.We do not provide CPU time reports for these six methods in the case of small n because they are less than one second of CPU time.
Below, we delve into the statistical significance of distinctions among the three heuristic approaches and the three GA algorithms.Initially, we applied an analysis of variance (ANOVA) technique utilizing a linear model on the AEP data for small 'n' or RPD data for large 'n' within SAS 9.4.The ANOVA, shown in Table 6 and 7, revealed that the model was significant for both small 'n' and large 'n'.Additionally, the algorithm and job-size (n) factors exhibited significance for both small 'n' with an F test value of 456.67 and five degrees of freedom and for RPD with an F test value of 50.13 and five degrees of freedom.It is also noticeable that the factors of weight (θ ) for the objective function and number of orders (m) are significant in Table 6 for small 'n', whereas they did not demonstrate significance for the large 'n' case shown in Table 7.The rationale behind this observation is that the effect of 'n' outweighs the values of weight (θ ) for the objective function or the number of orders (m) in the context of the large 'n' case.
The normality assumption of the linear model was validated for both the AEP and RPD data.This was evidenced by the Kolmogorov-Smirnov normality test, which produced p values exceeding 0.05 (with statistic values of 0.053336 for small 'n' and 0.058442 for large 'n').Additionally, the AEP and RPD Tukey grouping for means of algorithms were also tested.Figure 5 shows that the algorithms can be ranked into five groups.The five groups range from the worst to the best as DFOR, DPI, DBK, DFORGA, and (DPIGA, DBKGA) for the AEPs, shown in Fig. 5 (left) as well as the four groups from the worst to the best as The three heuristics (DPI, DFOR, and DBK) rely on predefined rules or strategies to find solutions quickly.While effective in some cases, they may get stuck in local optima and fail to explore the entire solution space.On the other hand, the three GAs (DPIGA, DFORGA, and DBKGA) employ population-based searches and a genetic-inspired selection, crossover, and mutation process, which allows them to explore a broader range of the solution space, potentially discovering better solutions.The advantages of genetic algorithms over heuristics in setup time optimization often stem from their ability to handle complex relationships, explore diverse solution spaces, and converge toward globally optimal solutions by using evolutionary mechanisms.Overall, it can be confirmed that the DBKGA is the best approach with the smallest value of AEP or RPD among all six proposed algorithms.

CONCLUSIONS AND SUGGESTIONS
In practical settings, a machine incurs a setup time whenever it switches from processing one product to another.Scheduling jobs while taking into account setup times and associated costs has attracted significant attention in both the manufacturing and service sectors, driving extensive research endeavors.While prior studies on customer order scheduling often focused on jobs to be processed across multiple machines, they frequently overlooked the critical factor of setup time.This study addresses a sequence-dependent bi-criterion scheduling problem, integrating considerations for order delivery.The primary aim is to pinpoint an optimal solution that minimizes the combined makespan and weighted sum of completion times for all specified orders.
To tackle this formidable challenge, we propose a branch-and-bound method incorporating two dominance properties and a lower bound.This approach is finely tuned to yield optimal solutions in scenarios involving a limited number of jobs, exemplified by cases with n =12, 16, 20, and 24.Demonstrating notable efficacy, this method excels for up to n =24 jobs, accommodating diverse combinations of order quantities and objective weights.For both small n (=12, 16, 20, and 24) and large n (= 60, 120, 180, and 240) tested instances, we introduce a tailored heuristic approach.This heuristic is complemented by three strategic techniques: pairwise interchange, extraction and forward-shifted reinsertion, and extraction and backward-shifted reinsertion (DPI, DFOR, and DBK).Additionally, we employ three heuristic-based genetic algorithms (DPIGA, DFORGA, and DBKGA) to obtain reliable approximate solutions.
Notably, it is evident that DBKGA performed significantly better than the five other proposed heuristics for both small and large test cases.A potential direction for future research could involve expanding the model to a sequence-dependent flow shop setting that incorporates order deliveries.Furthermore, it may be worthwhile to investigate extending this framework to minimize the total tardiness of all fulfilled orders.As another future topic, we can also consider that the machine can break down or extend this model to the multiple-machine setting.

Figure 1
Figure 1 Exploring the values of parameters for three GAs.Full-size DOI: 10.7717/peerjcs.1763/fig-1 . 3 present an evaluation of the performance of the heuristics and genetic algorithms in the context of larger n values.The data in Table Xu et al. (2024), PeerJ Comput.Sci., DOI 10.7717/peerj-cs.176313/23

Figure 5
Figure 5 The AEP and RPD Tukey grouping for means of algorithms (Alpha = 0.05).Means covered by the same bar are not significantly different.Full-size DOI: 10.7717/peerjcs.1763/fig-5 the last assigned job of O p , and J v is the last assigned job of O q , d Lu + +w q d vu and d Lu + d uv + max x∈π c {d vx } < d Lv + d vu + min x∈π c {d ux }, then σ dominates σ .Proof: In following, one will claim that (i) C max σ − C max (σ ) > 0, and (ii) The last inequality follows from the given condition d Lu + d uv + max x∈π c {d vx } < d Lv + d vu + min x∈π c {d ux }. (ii) For the orders O p and O q , [w q